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1. Introduction 

Lanchester-type* differential-equation combat models are an important 
tool for analyzing many important problems of military operations research. 
In such a combat model, a so-called attrition-rate coefficient represents 
the fire effectiveness of a particular weapon-system type against a particular 
target type, i.e. the weapon-system type's effective firepower against such 
a target. Time-dependent attrition-rate coefficients are used to model 
temporal variations in firepower on the battlefield. Thus, we see that 
time-dependent attrition-rate coefficients are important (and, in fact, 
essential [4-6]) for the quantitative analysis of hypothetical combat. 

Militarily realistic computer-based Lanchester-type models of quite 
complex military syst2ms have been developed for almost the entire spectrum 
of combat operations, from combat between battalion-sized units to theater- 
level operations. Nevertheless, a simple combat model may yield a clearer 
understanding of significant interrelationships that are difficult to per- 
ceive in a more complex model, and such insights can subsequently provide 
valuable guidance for more detailed computerized investigations. In this 
report we consider such a simplified variable-coefficient Lanchester-type 
model of combat between two homogeneous forces. 

For this variable-coefficient Lanchester-type model of combat 
between two homogeneous forces, different functional forms for the attrition- 
rate coefficiente lead to different mathematical functions being involved 
in representing and computing the force-level trajectories. In a previous 
paper [5] we have discussed the plausibility of the hypothesis that except 


for the special case of a constant ratio of attrition-rate coefficients, 


"so-called after pioneering work of F. W. Lanchester [3]. 


the solutions to such differential equations cannot be represented in terms 
of “elementary" functions of analysis. Thus, new transcendental functions 
arise in the study of combat modelled with time-dependent attrition-rate 
coefficients. In particular, we have previously introduced [5-6] so-called 
Lanchester-Clifford-Schlafli (LCS) functions for analyzing combat modelled 

with power attrition-rate coefficients with "no offset" (see Section 3 below). 


In the Appendix to this report is contained the most extensive 


set of tables currently available for the LCS functions: it contains tables 
of five-decimal-place values of the hyperbolic-like LCS functions F, (x), 

Hy 0%)» and Tox) (see Section 4 below) for 25 fractional values of 4 

(see Section 6 below). The main body of this report provides the theoretical 
and modelling background for the use of these tables. In particular, we 
examine a model of a constant-speed attack on a static defensive position 

and show how associated range-dependent kill rates give rise to time-dependent 
attrition-rate coefficients with "no offset." Numerical computations are 
presented to illustrate the use of the LCS functions for analyzing such 
“aimed-fire" combat. As a consequence of the availability of these tables, 


one can now study this variable-coefficient combat model almost as easily and 


thoroughly as Lanchester's classic constant-coefficient model. 


2. Variable-~Coefficient Lanchester-Type Equations of Modern Warfare. 


We consider combat between two homogeneous forces modelled by the 


following variable-coefficient Lanchester~type [3] (see [4,5]) equations of 


modern warfare 


; 
' 
’ 
H 


dx 
3 * -a(t)y with x(0) Xo > 

(2.1) 
ra = -b(t)x with y(0) = Yo? 


where t = 0 denotes the time at which the battle begins, x(t) and y(t) 
denote the numbers of X and Y at time t, and a(t) and b(t) denote 
time-dependent Lanchester attrition-rate coefficients, which represent the 
effectiveness of each side's fire. These coefficients depend on variables 
such as force separation, tactical posture of targets, rate of target acqui- 
sition, firing rate, etc. (see [4-7] for further details). Variable attrition- 
rate coefficients are used to model temporal variations in firepower on the 
battlefield. In any analysis of combat, moreover, we should use the above 
equations (2.1) only for x and y >O and, for example, set dx/dt = 0 
when x = 0, since negative force levels have no physical meaning. 
Mathematically, we assume that the attrition-rate coefficients a(t) 


and b(t) are defined, positive, and continuous for to <t<+o with 


ty < 0. We also assume that a(t) and b(t) € L(t),T) for any finite 
T> to: We further take a(t) and b(t) to be given in the form 
a(t) = k g(t) ; and b(t) = k h(t) ’ (2.2) 


where k, and ky are positive constants chosen so that a(t)/b(t) = ki/ 
when g(t) = h(t). We introduce the combat-intensity parameter Ay and 


the relative-fire-effectiveness parameter Ap defined by 


rs = vk ky, ; and dp =- k/k, n (2.3) 


From our assumptions about a(t) and b(t), it follows that, for example, 


T 
a(t) £ L(t)>T) implies f a(t)dt = +o. 
a) 


The X force level as a function of time may be represented 


as [5,6] 
x(t) = x, {C, (O)C, (t) a S,(0)s, (t)} - Yo" {c, (0)S, (t) - $,(0)C,(t)}, (2.4) 


where the hyperbolic~like general Lanchester functions (GLF) Cy (t) and 8, (t) 


are linearly-independent solutions to the X force-level equation 


2 
dx 1 da dx 
si - tats al ae a(t)b(t)x = 0, (2.5) 
with initial conditions 
Cy (t) #=l, $y (to) = 0, 
(2.6) 
{1/a(to)} dC /dt(t,) = 0, {1/a(t,)} 48,/dt(t,) = 1/va, f 


Here to denotes the largest finite time at which a(t) or b(t) ceases 
to be defined, positive, or continuous. The Y force level as a function 
of time is given by a similar expression, with cy (t) and Sy (t) being 
analogously defined for the corresponding Y force-level equation. 


It is sometimes convenient to introduce the new independent variable 


t defined by 


Pace Paes ee 


t 
t= { Va(s)b(s) ds . (2.7) 
t , 
0 


It is readily seen that the transformation t= t(t) is well defined and 


invertible. Let us denote 1(0) as To: We observe that to <0 implies 
that 1) > 0. If we denote the “average intensity of combat" as 
Yva(t)b(t) , then 
ee eter oe ae t 
Ya(t)b(t) t = {(1/t) { Ya(s)b(s) ds}t =1- Tj)? (2.8) 
0 
The substitution (2.7) transforms (2.5) into 
dx - ey {2 Q R(x} as .. x = 0 (2.9) 
a. @wie dt , : 
dt 
with initial conditions 
x(t) = Xp > and {1/¥R(x5)} dx/dt(t)) = -¥q » 


where R(t) = a(t)/b(t). 


3. Combat Modelled with Power Attrition-Rate Coefficients. 
The above equations (2.1) basically apply to "aimed-fire" combat 
when target-acquisition times do not depend on the numbers of targets available 


(see [5,6] for further details). A large class of tactical situations of 


coefficients [5-7] 


i 
interest can be modelled with the following general power attrition-rate | 


a(t) = k(t + c)", and b(t) = k(t +C+ 4°. (3.1) 


where A and C > 0. We will call A _ the offset parameter, since it allows 
us to model (with j and v > 0) battles between opposing weapon systems 
with different maximum effective ranges (see [5,6]). We will call C the 
Starting parameter, since it allows us to model (again, with wu and vy > 0) 
battles that begin within the maximum effective ranges of the two opvosing 
systems. We observe that for the general power attrition-rate coefficients 


(3.1) we have = -C, and u and vy must be > -1l in order that a(t) 


*o 
and b(t) € L(t),T). 
The above nomenclature is motivated and possible applications of our 


work are indicated by considering S. Bonder's model of the constant-speed 


attack on a static defensive position (see [4-7] for further details) 


dx d 

= -a(r)y, and = = -B(r)x , (3.2) 
where r denotes the range between opposing forces, and a(r) and 8(r) 
denote range-dependent attrition-rate coefficients. Range is related to time 
by 


r(t) = R, -vt, (3.3) 


0 


where Ry denotes the opening range of battle and v > 0 denotes the 
constant attack speed. For example, let us consider the constant-speed attack 
of a homogeneous Y force against the static defensive position of a homo- 


geneous X force. Figure 1 diagrammatically portrays this situation. 
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The basic idea is that force separation, i.e. range between the opposing 
forces, changes over time, and the fire effectiveness of, for example, 
a single Y firer, denoted as a(r), depends on this force separation. 

In many cases of tactical interest, we may model the fire effective- 


ness of, for example, the Y weapon system (as a function of range) with 


r \" 
% l- R, for O<r<R, 
a(r) = (3.4) 
0 for R <r, 
se 


where Ro denotes the maximum effective range of the Y weapon system and 
u > 0, Here yw is used to model the range dependency of Y's attrition-rate 
coefficient (see Figure 2). We model &8&(r) similarly, with corresponding 
quantities R, and v_ being analogous to R. and wu above. 


If we use (3.3) to eliminate range r from (3.4), we obtain 


dx 
— 
(3.5) 


ao 
at b(t)x . 


where the time-dependent attrition-rate coefficients a(t) and b(t) are 


given by (3.1). It follows that the offset and starting parameters are given 


by 
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R, - R R -R 
rp oe a ee ee 6.6 


and that 


Vv u Vv Ww 
k, = 29 (3) F and k, = By (x) : (3.7) 
a B 


vc observe that A and C >0 if and only if R, > R. > Ro: By consider- 
ing (3.6) and Figure 3, the reader should have no trouble in understanding 
our terminology for A and C. 

When the offset parameter is equal to zero (i.e. A = 0), then the 


coefficiente (3.1) reduce to 


a(t) = k(e+c)", and b(t) = k(t +c)” ’ (3.8) 


We will refer to (3.8) as power attrition-rate coefficients with "no offset." 


As we have seen above in Bonder’s constant-speed attack model, these 
coefficients model, for example, combat between weapon systems with the same 
maximum effective range so that there is no "offset" in the "reaching out" 
of the weapon systems against each other in combat (again, see Figure 3). 
For these coefficients (3.8), the transformed X force-level equation (2.9) 


becomes 


2 
oxy @al & yo, (3.9) 
de t ° dt 


with initial conditions 
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q utv+t2] ” . 
2d 
I (u+v+2)/2 
ce(peta) ero ’ (3.11) 
and 


2r 
‘o> ( I pa aa) 


utv+t2 


Let us observe that 0 <q <1 when y» and v >-1l. Furthermore, q > 1/2 
if and only if dR/dt <0, i.e. R(t) is a strictly decreasing function of 


time. 


4. Lanchester-Clifford-Schlafli (LCS) Functions. 
Consider the function F,(x) defined by the power series 


(= ra) J x/2)"* 
F (x) = [(a aes : (4.1) 
a k=0 k! T(k+a 
For a #0, -1, -2,... the radius of convergence for F,(x) is infinite 


by the ratio test for convergence of power series (2). Hence, F (2) is 


an entire function of the complex variable z= x + iy, with an essential 


12 


a 


singularity at the point at infinity. Now consider the function H (x) 


defined by the infinite series 


2(k+a) 


. , x/2) 
H (x) l(a) oo STS CEarEESD ‘ (4.2) 


Observing that 
H (x) = (1/a) (x/2)2% Fig) > (4.3) 


we see that for a>0O the infinite series (4.2) is uniformly convergent 
on compact subsets of the complex plane. From (4.3) one can readily deduce 


the recursive relation 


2 
+ |B F 


F(x) = F CK: (x) . (4.4) 


ot 


We will call the functions F (x) and H, (x) Lanchester-Clifford-Schlafli 


(LCS) functions (see Note 10 on pp. 66-67 of [5]). Other properties are 
readily deduced and are given in Table I. 
The function F(x) satisfies the linear second-order ordinary 


differential equation 


dF 
a 2a-1 a 
—— + ae )}—- Fr. #0, (4.5) 


with initial conditions 


Table I. Properties of the LCS Functions F (x) and Hy (x). 


1-2 
lL. dP /dx = (x/2) 7H, (x) 
2. di /dx = (x/2)7*"*F, (x) 


3. F(x) Fy, (*) ~ Hy (x) Hy _, Gx) -=1 #£Vx 


where a is not an integer (including zero) 


4. F, (x=0) =1 


5. H, (x=0) = 0 for a>0Q 


6. dF /dx(x=0) = 0 


1-20 


7. {(x/2) dH /dx} 29 =1 


8. Fij2® = cosh x 


9. By 72‘) = sinh x 


dF 
z - oe oe 
F (0) 1, and ax (0) =0, 
while Hy Gx) satisfies 
2 
di dH 
a 2a-1 a 
fey ae oe ee ci 
dx 


with initial conditions 


1-2 dH 
H(0)=0, 3 and >) aH, |  e 
is  x=0 


Thus, {F, Hy.) is a fundamental system of solutions to 


2 
—- + BE = F=0, (4.7) 


with Wronskian WF. H, _,.) = (x/2)1729, It follows that the GLF for the 
X and Y force-level equations for combat modelled with the attrition-rate 
coefficients (3.8) are given by 


Ay 


Shetd Ho (x(t) . (4.8) 


C, (t) = Fi (r(e)), 8, (t) -( 


Ay -- 


Cy (t) = Fy (r(t)), Sy(t) (sts Hy (x(t) ‘ (4.9) 


where p= 1-q. If we define 


T(x) = Hy _, (0) /F, (x) ’ (4.10) 
then 
S(t) Ay ‘oe H (1t(t)) 
KO "Ca “\iev +2 PCD ” ead? 
or 
( Ay - 
Ty (t) ek FF es Ty Cr(e)) > (4.12) 


where Ty (t) denotes a hyperbolic-like GLF, which corresponds to the 


hyperbolic tangent. Observing that for wu, v > -l, lim t(t) =+°, 


TH+ oO 


we see that T(x) is a strictly increasing function of x on the interval 


[O, + ~) and 
0 < T(x) < C=) for O<x<+-, (4.13) 
with 
. -G-a) 
bon T, (x) T(a)’ (4.14) 


since by the results of Taylor and Comstock [7] the parity-condition parameter 


Q* = Q*(u, v, C = 0) is given by 


x 2q-1 
engender a I r(p) 
ron 5 Q*(u, v, 0) (; +v+t :) T(q) ° (4.15) 


We recall that Taylor and Comstock [7] have introduced the so- 
called parity-condition parameter Q* as the value (or range of such values) 


for the initial condition Q to the initial-value problem 


meme aro 


\ 

t | 
pa | 
ee | 


tain a(t)E, with E, (to) _ ; 

R (4.16) 
dE Es as 
, ade Yrp b(t)EY with E(t.) =Q, 


such that EY (t3Q*) and By (t3Q*) > 0 for all t > to: In other words, 

Q* is the value of Q in (4.16) above such that neither EY nor Ey 

ever become zero. In this case, both Ey (t3Q*) and Ey (t3Q*) are positive, 
strictly decreasing functions, similar to decreasing exponentials. Thus, 

we may call Q* "the Y equivalent of an X force of unit strength," 

since the forces are "at parity," with neither force being annihilated in 
finite time. Taylor and Comstock have shown that for either a(t) ¢ L(0,+ ~) 


or b(t) ¢ L(O0, + =), then Q* is unique and given by 


S(t) 
Pocih GC) * P (4.17) 


The significance of the parity-condition parameter Q* is that it allows 


us to predict force annihilation as the following theorem shows. 


THEOREM 1 (Taylor and Comstock [7]): Assume that either 
a(t) # L(O, +=) or b(t) ¢ L(O, +”). Then the xX force 


will be annihilated in finite time if and only if 


C,(0) - Q*S, (0) 
Vip HCO) = 5, OY ‘ (4.18) 
17 


5. Use of LCS Functions for Analyzing Combat. 
The Lanchester-Clifford-Schlafli (LCS) functions F(x) and H, (x) 


are useful for analyzing "aimed-fire" combat (see Section 3 above) modelled 
with the power attrition-rate coefficients with "no offset" (3.8), which 


we rewrite here as 
a(t) = k(t +c)", and b(t) = k(t +c)”. (5.1) 


In other words, the LCS functions arise in solving the differential combat 
model (2.1) with attrition-rate coefficients (5.1). In order that both 
a(t) and b(t) € L(ty,T), we must have uw and v>-l. Military situations 
modelled by these equations have been discussed in Section 3 above, e.g. 
combat between two weapon systems with the same maximum effective range. 
For such combat, the LCS functions may be used to 

(1) compute force-level declines, 

(2) predict force annihilation, 
and (3) predict the time of force annihilation. 

Let us now see how the LCS functions may be used to obtain the 
above information about force-level declines and force-annihilation pre- 
diction. According to (2.4), (4.8), and (4.9) above, the X force level 


is given by 


x(t) = Xo fF (to) F(t (e)) - H(t) H(r(e))} 


Ay 2q-1 
- 99% \oeS (Fy (to)H (Ce) -H(rg)F (rCe))}, (5.2) 


18 fee | 


where q is given by (3.10), p = l-q, and ;(t) is given by (3.11), which 


we rewrite as 


2 
I (u4v+2) /2 
T(t) = ( sete) (t + C) , (5.3) 


The time to annihilate the xX force* is determined by x(t) = 0, and 


it follows that 


rr, A, \a-p 
gv. Se p'o t%o alae ey 3% 
a Ore 7* 


x =P ’ (5.4) 
XH (to) + Yq Mp =) Fy? 
where from (4.10) 
¥, er€ed) = Hy (1 (t))/F (r(t)) ’ (5.5) 


and we recall that p+q#1. It follows that the time to annihilate X, 


r* 


a’ is given by 


"Te we multiply the first equation of (2.1) by y, the second by x, add, 
and integrate, we obtain 


t 
x(t) y(t) = XoVo ~ J {a(s) y*(s) + b(s) x7 (8) }ds ; 
0 


which shows that x(t) and y(t) can have at most one finite zero, Hence, 
if x(tX) = 0, then we know that y(t) > 0 for all t > 0. 


(1,) / ( “1 ib (1) 
a) i | XP (te) tp Mp | —aaS H(t 
Xe cl) g-2 | “orp to o’*R \ptv+2 p-‘0 2 oY os gs 


( x i 
a 
XH t9) . Yo YR utv+t2 Fy) 


Taylor and Comstock [7] have shown that Tf is strictly increasing and 


satisfies (see also (4.12) above) 
7 < T(p)/T(q) , (5.7) 


where p = l-q. It follows that in order for X _ to be annihilated in finite 


time, the right-hand side of (5.4) must be less than [(p)/I'(q). Let us 


observe that for ty = -C = 0, (5.4) simplifies to 
P-q 
T(x (e)) - = eet) ‘ (5.8) 
Yo "*R 


Thus, we have proved the following theorem concerning force~ 


annihilation prediction. 


THEOREM 2: Consider combat between two homogeneous forces 
modelled by (2.1) with power attrition-rate coefficients 
(5.1). Assume that yp and v > -1 and that the above 
equations hold for all time. Then the X _ force will be 


annihilated in finite time if and only if 


| ee ee 
PCa) | XoF, (To) * Yo VAR acted H(t) 


q-P 
| % ) 
< TP) ) XH (tq) + ¥9 MQ EV Ee Fylto) ( ° (5.9) 


where q = (v+1)/(u +v +2) and p=1-q. For to = 0 


(i.e. C= 0 so that 1, = 0), X will be annihilated in finite 


0 
time if and only if 
Xo r(p) dy ae 
Yo : T'(q) re utv+t2 : (5.10) 


6. Tabulation of LCS Functions. 
This report contains the most extensive set of tables of the 
Lanchester-Clifford-Schlafli functions currently available. The Appendix 


contains tables of five-decimal-place values of the hyperbolic-like LCS 


functions F, (x), Hy _,(), and TG) for various values of the argument x, 
namely x = 0.00 (0.01) 2.00 (0.1) 10.0, and a = 1/2, 1/3, 2/3, 1/4, 3/4, 
1/5, 2/5, 3/5, 4/5, 2/7, 3/7, 4/7, 5/7, 4/9, 5/9, 3/11, 5/11, 6/11, 8/11, 
5/13, 8/13, 5/17, 12/17, 5/21, and 16/21. These values of the index a 
correspond to u, v = 0, 1/4, 1/2, 1, 1 - 2, and 3 in (3.8) and allow one 
to analyze, for example, a fairly wide variety of range capabilities for 


weapon systems in the constant-speed-attack model of Section 3. These 


| 


4 F <e 
ees 


tables have been calculated by the recursive means given in Section 8 of 
[5]. A less extensive tabulation (namely, for a = 1/2, 1/3, 2/3, 1/4, 3/4, 
1/5, 2/5, 3/5, 4/5, 3/7, and 4/7 corresponding to wu, v= 0, 1, 2, 3) is 
to be found in a companion report [8]. 

A representative tabulation of the hyperbolic-like LCS functions 
FL, Hy)» and Tx) is given in, for example, Tabies 8A and 8B of 
the Appendix for a = 3/5. The values of the argument x are the same as those 
used for the tabulation of the hyperbolic functions by Abramowitz and Stegun 
[1]. We observe from Table 8B and (4.13) that the limiting value of 
T(x) as x*+© (here a = 3/5) is quickly reached, with three-decimal- 
place accuracy already attained for x = 4.5. Moreover, the use of these 
tables (specifically, Tables 8A and 8B of the Appendix) for combat analysis 


is illustrated in the next section. 


7. Numerical Examples 


In this section we examine a couple of numerical examples to show 
some of the insights that may be gained into the dynamics of combat between 
two homogeneous forces from our results (see also [6]). These examples 
illustrate the use of the LCS functions Fix), Hy), and 7) for 
analyzing "aimed-fire" combat modelled with the power attrition-rate coefficients 
with "no offset" (5.1). As in [4-7], we consider S. Bonder's model (3.2) 
for the constant-speed attack against a static defensive position. We will 
focus on the use of the LCS functions for predicting force annihilation, 
since the computing of force-level trajectories with Lanchester functions 


is adequately handled elsewhere (see [4-5]). ae 
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Let us accordingly consider the constant-speed attack of a homogeneous 
Y force against the static defensive position of a homogeneous xX force 
(see Section 3 above for further modelling details, especially Figure 1). 
For our numerical computations, we assume that the fire effectiveness of the 


Y weapon system varies linearly with range, i.e. 


a) (1 - =} for O<r<R, 


a(r) = (7.1) 


0 for R_ <r, 
os 


and that the fire effectiveness of the X weapon system varies quadratically 


with range, i.e. 


(1 - x) for O<r<R 


B(r) = (7.2) 


with R. = Rp» i.e. both weapon systems have the same maximum effective range. 
In other words, 1, = 1 in (3.4) and v= 2 for 8(r). We consider a battle 
modelled by the input data given in Table II. In terms of time as the 
independent variable, the attrition-rate coefficients (7.1) and (7.2) become 


via (3.3) 


a(t) = k(t +€) and b(t) = k(t + c)* , (7.3) 
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Table II. Input Data for Numerical Examples 


wel, ve=2 


= 0.06 X casualties/minute/Y firer 


= 0.6 Y casualties/minute/X firer 


R = Re = 2000 meters 


v = 5 miles/hour 


Av re 2 
Cc: ——— Fy k, = s, : and k, = t0( zt) « (7.4) 


From the input data given in Table II, we compute the parameter values shown 

in Table III, since the transformed X force-level equation is given by (3.9) 
with q= (v+1)/(u +v+ 2), p = leq, » = 1, and v= 2. Thus, the X force 
level may be computed with Fi) and Hy) with a= q = 3/5. Force- 
annihilation prediction involves the limiting value of T(t) = Hy _ (D/P) 
as t ++, From Table 8B of the Appendix and Table III, we note the predicted 
agreement between [(1-«)/[(a) and the limiting value of 1, G) “as «** * 
[recall (4.13)] for a= q = 3/5. We now consider two cases: (I) Ry = 2000 
meters, and (II) R, = 1250 meters. 


0 


When Ry = 2000 meters (see Figure 3 of [4]), we have C =O and 


tT * 0. The maximum time that the battle can last is os Ry/v = 14.91 minutes, 
since at this time the attackers reach their final objective, i.e. the 

defender's position (again, see Figure 1). We now consider the qualitative 
behavior of the » = 1, v= 2 force-level trajectory shown in Figure 3 of 


[4]. Theorem 2 tells us that the X force can be annihilated in finite time 


if and only if 


Xo r(p) dy q-P 
Yo * F(a) OR (; + ry) , (7.3) 


where q = 3/5 and p= 1-q. Using the numerical values in Table III, we 


compute from (7.3) that the X force can be annihilated in finite time if 


and only if 


Table III. Parameter Values for Numerical Examples 


k, = 4.0233 x 16" X castaltios/ininute)" © fiver 


kK = 2.6979 1073 Y casualties/ (minute) ”/x firer 


x 


p = 2/5, q = 3/5 


l'(p)/T(q) = 1.48951 


A=0 


— < 0.420. (7.4) 


When the X force can be annihilated, its annihilation time is given by 


(5.8), which we rewrite here as 


x Xo Ay P-4 

T (t(t_)) = (sz) ° (7.5) 

q a utv+2 
Yo Vie 

where 
2A 
I (utvt2) /2 : 
t(t) = ie t 4 (7.6) 


Thus, for the numerical values given in Table III, the time of annihilation 


of the X force is given by 


x Xo 
T (r(tQ)) = 3.544 %” (7.7) 


with q = 3/5. We will now illustrate further computations for x, = 10 and 


0 
Yo * 30. From (7.4) we see that the X force can be annihilated in finite 
time (but we must verify that e < tax)’ In this case (7.7) becomes 
Taxes) = 1.18122 . (7.8) 


We must now determine r(t%) such that r(t%) - 1) (1.18122) by using inter- 


polation methods and Tables 8A and 8B. From Table 8A, we find 
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ae om ee 


To = 1.18172 for = 1.01 


¥,() = 1.17630 for = 1.00 
so that using linear interpolation, we obtain 
1(t2) = 1.009 , (7.9) 
whence use of (7.6) yields 
ty = 14.24 minutes, (7.10) 


which is less than a = 14.91 minutes so that the defending X force is 

indeed annihilated before the attacking Y force reaches its final objective. 
Since r(t) = Ry 
annihilation of the X force is 


- vt, we find that force separation at the instant of 


r = 89.8 meters . (7.11) 


Further results may be computed in a similar fashion and are given in 


Table IV. 

When Ry = 1250 meters (see Figure 3 of [5]), we have C = 5.5923 
minutes, tT)" 0.0975, and as Ro/v = 9,32 minutes. In this case 
Theorem 2 tells us that the X force can be annihilated in finite time 


1f and only if 
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Table IV. Annihilation of the X Force as a Function 


of the Initial Force Ratio for Ry = 2000 meters 


(x9/¥q) a (minutes) ¥ (meters) 
0.333 14.24 89.8 
0.250 11.61 443.2 


0.200 10.19 633.2 


Pq) 
(t)) -F 


A - H 
% (see) 38 ‘i 9} avis 
y R\utv+2 T'(q) r . P 
0 ig) - Fs H(t) } 


with q = 3/5 and p = l-q. Using linear interpolation, we obtain from 


Tables 7A and 8A of the Appendix that for the numerical values of Table III 


FM) = 1.006 , H(t») = 0.044 , 
(7.13) 
Ff) = 1.004 , H(t) = 0.223 , 
so that (7.12) says that the X force can be annihilated if and only if 
*o 
— < 0.382 . (7.14) 
Yo 
When the X force can be annihilated, its annihilation time is given by 
(5.4), which we rewrite here as 
P-q 
x r 
0 I ) 
a rae F (t,) + H(t 
, a (=. p 0 io 
= 7.15 
Ty (ree) Xp , P-q ’ ( ) 
FG) + — ah.) H(t) 
q 0 utv+2 q' 0 
Yo" AR 


whence for the data of Table III 


x. 3365, + 0.223 
T,(t(e.)) = 0.156u, + 1. , (7.16) a 


where Uy * Xo/Yo- Let us also record here that (3.11) yields 


t= (test zie See (7.17) 


2/ (ut+v+2) 
ame. 
We will again illustrate further computations for X * 10 and Yo * 30. 
From (7.14) we see that the X force can be annihilated in finite time (but 
x 


again we must investigate whether or not = Sie? In this case (7.16) 


becomes 
xX 
Fy (re.)) = 1.33651 , (7.18) 
whence Table 8A of the Appendix and linear interpolation yield 
x 
1(t,) = 1.397 , (7.19) 
so that by (7.17) 


2 = 10.63 minutes . (7.20) 


a? We see that the attacking Y force 


Since _.* R/V = 9.32 minutes < t 
overruns the defender's position before annihilation of the X force occurs. 
Thus, the battle ends with x," x(t,) > 0 and Ye ” 0 at t= iS 
9.32 minutes. Corresponding to te = 9.32 minutes is t, *™ 1.1318, and 


then Table 8A of the Appendix yields 
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F(t = 1.1318) = 1.589 , M Cter326) = 1.973 , (7.21) 


whence via (2.4), (4.8), (4.9), and (7.13) we obtain 


X, = x(t ¢) =x(r = 0) = 1.35. (7.22) 


Some further numerical results are given in Table V. Again, these parametric 
results should be contrasted with the single wu =1, v= 2 force-level 


trajectory shown in Figure 3 of [5]. 


8. Final Remarks 

In the previous section above, we have seen how the LCS functions 
allow one to conveniently obtain much valuable information about the model 
(2.1) with power attrition-rate coefficients (3.8) without having to explicitly 
compute the entire force-level trajectories. Previously we were limited to 
computing only force-level trajectories (see [4-5]). With the availability 
of these tabulations of LCS functions (see the Appendix of this report and 
[8]), we can now tell who is going to be annihilated and when this event will 
happen without having to compute the trajectories. Not only did we answer 
questions about the qualitative behavior of the model (e.g. force annihilation) 
for specific values of, for example, initial force levels but also for a 
range of values of the initial force ratio (i.e. parametric analysis of 


model behavior). - 
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Table V. Annihilation of the X Force as a Function 


of the Initial Force Ratio for R, = 1250 meters 


0 
(X9/¥o) r* (minutes) rx (meters) 
a a 
0.333 10.63 +? 
0.250 7.56 235.9 
0.200 6.17 422.8 


_ ae 9.32 minutes and x. * x(r=0) = 1.35. 


The results of this report may be used for other parametric 
analyses, e.g. parametric dependence of battle outcome on attrition-rate 
coefficients. Thus, the contents of this report allow one to develop 
important insights into the dynamics of combat between two homogeneous 
forces with temporal variations in fire effectiveness. With the avail- 
ability of tabulations of the LCS functions, one can now analyze such 
combat modelled by the power attrition-rate coefficients (3.8) with some- 
what the same facility as he can for the constant-coefficient case and thus 
aid in parametric analyses. For further discussions of the significance 
of such results for military operations research, the reader is directed 


to [6] and [7]. 
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APPENDIX: Tabulation of the LCS Functions F(x), H,_,@), and T(x) for 
a = 1/2, 1/3, 2/3, 1/4, 3/4, 1/5, 2/5, 3/5, 4/5, 2/7, 3/7, 
4/7, 5/7, 4/9, 5/9, 3/11, 5/11, 6/11, 8/11, 5/13, 8/13, 


5/17, 12/17, 5/21, and 16/21. 
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